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The scheme of entanglement-assisted quantum error-correcting (EAQEC) codes assumes that the 
ebits of the receiver are error-free. In practical situations, errors on these ebits are unavoidable, 
which diminishes the error-correcting ability of these codes. We consider two different versions of this 
problem. We first show that any (nondegenerate) standard stabilizer code can be transformed into 
an EAQEC code that can correct errors on the qubits of both sender and receiver. These EAQEC 
codes are equivalent to standard stabilizer codes, and hence the decoding techniques of standard 
stabilizer codes can be applied. Several EAQEC codes of this type are found to be optimal. In 
a second scheme, the receiver uses a standard stabilizer code to protect the ebits, which we call a 
"combination code." The performances of different quantum codes are compared in terms of the 
channel fidelity over the depolarizing channel. We give a formula for the channel fidelity over the 
depolarizing channel (or any Pauli error channel) , and show that it can be efficiently approximated by 
a Monte Carlo calculation. Finally, we discuss the tradeoff between performing extra entanglement 
distillation and applying an EAQEC code with imperfect ebits. 



I. INTRODUCTION 

The theory of quantum error correction is important 
for both quantum computation and quantum communi- 
cation Quantum stabilizer codes are the most ex- 
tensively studied quantum codes [7, 8 , and have the ad- 
vantage that their properties can be analyzed using group 
algebra. Quantum stabilizer codes are closely related to 
classical linear codes, and can be obtained by the CRSS 
and CSS code constructions from weakly self-dual clas- 
sicalcodes [7HT2]. When entanglement between sender 
and receiver is available, a new error correction scheme 
becomes possible: entanglement-assisted quantum error- 
correction. This coding scheme (EAQEC codes) has the 
advantage that it allows any classical linear code, not nec- 
essarily weakly self-dual, to be transformed into a quan- 
tum code [13]. In addition, EAQEC codes can increase 
both the transmission rate and error-correcting ability 
O Also, some problems or limitations in quan- 
tum LDPC codes and turbo codes can be solved using 
EAQEC codes [l5l[16]. 

In EAQEC codes [13], it is assumed that the sender 
(Alice) and the receiver (Bob) share some pairs of qubits 
in maximally-entangled states (ebits) before communica- 
tion, and the quits on Bob's side are subject to no error. 
The quantum codes are designed to cope with the noisy 
channel Ma that Alice uses to communicate with Bob. 
The properties of EAQEC codes in this case are studied 
in [131 HH [17] . However, noise (such as storage errors) 
can occur on Bob's ebits in practical situations, which 
is believed to degrade the performance of the quantum 
codes. 
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Assume the errors occurring on Bob's qubits are de- 
scribed by a noise process Mb- Wilde and Hsieh 
addressed this question with a channel-state coding 
protocol in quantum Shannon theory and determined 
the channel capacity when entanglement is not perfect 
[TS] , They also performed simulations of entanglement- 
assisted quantum turbo codes with the depolarizing chan- 
nel when Bob's ebits also suffer errors [15]. Wilde and 
Fattal simulated the performance of an entanglement- 
assisted Steane code for fault tolerance fT9|. 

In this article, we discuss two coding schemes to han- 
dle the problem when the ebits of Bob are not per- 
fect. Shaw et al. described a six-qubit EAQEC code 
with one ebit that is equivalent to Steane's seven-qubit 
code, and can correct a single error on either Alice's 
or Bob's qubits [20 . The entanglement-assisted Steane 
code, constructed by Wilde and Fattal, is also equiva- 
lent to Steane's seven-qubit code [19j. Similarly, Bowen's 
entanglement-assisted code [2T1 is equivalent to the five- 
qubit code [5] [22] and can correct an error on one of 
Bob's qubits. These three examples motivate the follow- 
ing idea: there are EAQEC codes that are equivalent to 
standard stabilizer codes, and hence can correct errors 
on both Alice's and Bob's sides. We show how to obtain 
an EAQEC code from a (nondegenerate) stabilizer code. 
Several EAQEC codes from this scheme are found to be 
optimal. We say a quantum code is optimal if the mini- 
mum distance of this code achieves an upper bound for 
fixed numbers of information qubits and physical qubits. 
These EAQEC codes will have better performance than 
their equivalent stabilizer codes when the storage error 
rate is less than the channel error rate. 

In the second scheme, Alice uses an EAQEC code to 
encode her information qubits and Bob uses a standard 
stabilizer code to protect his halves of the ebits. The 
combination of an EAQEC code and a stabilizer code 



2 



is called a combination code, and it can be treated ei- 
ther as a single stabilizer code, or by using two sequen- 
tial decoders. EAQEC codes that are not equivalent 
to standard stabilizer codes generally have higher error- 
correcting ability on Alice's qubits and are suitable for 
this scheme. 

Minimum distance of a stabilizer code is used as a mea- 
sure of how good a code is without considering the de- 
tails of the noisy channel model. However, minimum 
distance might not always be the best measure, for a 
quantum code may be able to correct many error opera- 
tors of weight higher than that indicated by the minimum 
distance. In particular, there is no general definition of 
minimum distance for the variant coding schemes in this 
article. A perhaps more suitable merit function is the 
channel fidelity 23 IS], which compares the similarity 
of the modified quantum state with the original quan- 
tum state. However, the calculation of the channel fi- 
delity depends on the channel and has an exponentially 
increasing complexity. We derive a formula for the chan- 
nel fidelity of a quantum stabilizer code over the depo- 
larizing channel, which facilitates its computation. The 
channel fidelity also can be well approximated by a lower 
bound when the depolarizing rate is small. Furthermore, 
Monte Carlo methods can often efficiently approximate 
the channel fidelity [25 . 

Another natural question arises in EAQEC codes. The 
perfect entanglement shared between sender and receiver 
will in practice be generated from a process of entan- 
glement distillation [5, 26 or a breeding protocol [27l. 
It is known that entanglement distillation with one-way 
classical communication is equivalent to a quantum error- 
correcting code [5 . Since we can also communicate using 
an EAQEC code that is robust to imperfect ebits, we dis- 
cuss whether it is always necessary to do entanglement 
distillation before communication, and how much. 

This paper is structured as follows. Basics of stabilizer 
codes and EAQEC codes are reviewed in the next sec- 
tion. Basic criteria for an EAQEC code to be capable of 
correcting errors on Bob's qubits are analyzed in Section 
HI We discuss the first coding scheme with imperfect 
ebits in Section \P7\ and the second scheme in Section 
[V} The formula for the channel fidelity over the depolar- 
izing channel is derived in Section |VI[ and the channel 
fidelities for quantum stabilizer codes and EAQEC codes 
are given in Subsection |VI A| and |VIB[ respectively. We 
discuss Monte Carlo simulations for the channel fidelity 
in Subsection |VI C[ In Section [Vll[ we compare the per- 
formances of different coding schemes in terms of the 
channel fidelity. Entanglement distillation is discussed in 
Section IVIIII The conclusions follow in Section HXl 



II. BASICS 

Let H be the state space of a single qubit. Suppose Al- 
ice sends a /c-qubit state to Bob by using an [[n, /c, d]] 
quantum stabilizer code that encodes the /c— qubit infor- 



mation state in a 2^-dimensional subspace of the n- 
qubit state space H'^^^ fixed by a stabilizer group S with 
some minimum distance d. The stabilizer group S is an 
Abelian subgroup of the n-fold Pauli group Qn^ with n — k 
generators ^i, ^2, • ' ' 7 9n-k^ and does not contain the neg- 
ative identity operator — /. Suppose /J^; is an n-fold uni- 
tary Clifford encoder^ which leaves the n-fold Pauli group 
Qn invariant under conjugation. The encoded state is 
Ue (1^^)10)*^"^"^). For convenience, let the stabilizer gen- 
erators be gi = UEZi-^kU^j for z = 1, • • • , n — /c, where 
the subscript i -\- k of the Pauli operator means that the 
operator is on the {i + k)-th qubit. The error-correction 
condition for stabilizer codes says that {Ei} is a set of cor- 
rectable error operators in Qn if EjEj ^ X{S)\S [51112]. 
where JV{S) is the normalizer group of S in Qn- (Since 
the the overall phase of a quantum state is not impor- 
tant, we consider errors of the form Mi (g) • • • (8)M^, where 
Mj e {/, X, r, Z} for j = 1, • • • , n.) This implies the def- 
inition of the minimum distance d of a stabilizer code to 
be the minimum weight of any element in A/'(*S)\*S, where 
the weight wt(^) of g G Gn is the number of components 
of g that are not equal to the identity. 

After Bob receives the noisy quantum state, he does 
the following three steps to recover the information state 
syndrome measurement, correction, and decoding. 
He first applies a series of projective measurements with 
projectors 



n 



on the output state of the noisy channel. Here, s = 
S1S2 ■ ■ ■ Sn-k is a binary (n — /c)-tuple that represents the 
error syndrome. The error syndrome s of an error op- 
erator E has Sj = if E commutes with gj and sj = 1 
otherwise. We sometimes represent s as a number with 
binary form si ■ ■ ■ Sn-k- Note that the P^'s are orthog- 
onal to each other, and Pq is the projector on the code 
space, that is. 
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PoUe {\^) ® 10)®"-'=) = Ue {\^P) ® |0) 



for any /c-qubit state |?/^). Given a stabilizer group 
5, there are 2"^"^ distinct error syndromes. For each 
nonzero error syndrome s, we choose a Pauli operators Eg 
(not in tS), whose error syndrome is s; for the error syn- 
drome s = 0, we choose £^0 = E The error operators Es 
are called syndrome representatives. If the measurement 
result is 5, the correction operator Cs = Es is applied, fol- 
lowed by the decoding unitary operator U^j. Finally, Bob 
throws away the ancilla qubits, which is the same as ap- 
plying a partial trace over the ancilla qubits. We define 
a set T containing the syndrome representatives {Eg}. 
Then |T| = 2'^~^ and T is a set of correctable error oper- 
ators. Note that the choice of T determines the decoding 
process. {S determines the encoding process.) In fact, we 
have many more correctable error operators than T. For 
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g £ the operation of Egg and Eg on the encoded state 
are the same (^5^t/£;|^)|0)^^-^ = EsUE\i^)\0)^''-^) and 
can be corrected by the the same correction operator Cg . 
The error operator Egg is cahed a degenerate error of 
Eg. The set T X 5 = {hg : h e T, g e S} is 8i correctable 
set of error operators that satisfies the error correction 
condition. 

Lemma 1. For a given T, the error operators in Tx^S are 
correctable and are the only correctable error operators. 

Therefore, we have a total of 4^~^ correctable n-fold 
Pauli operators. Note that a code can have different and 
inequivalent sets T. 

In this article, error processes are modeled by the de- 
polarizing channel T^^ independently operating on n 
qubits, where 

Tp{p) = (1 - ^p)p + ^{XpX + YpY + ZpZ), 

with depolarizing rate p for < p < 1 and p the density 
operator of a single qubit. Since the operation elements 

of the depolarizing channel Tp are ~ \v^-> \l\p^-> 

'^F, and yJ\pZ^ the operation elements of T^^ are 

{^/PiEi}^ where Ei is a Pauli operator in the n-fold Pauli 
group Qn and Pi is the probability that error Ei happens. 
If Ei is of weight then pi is 
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Now assume c maximally-entangled states 1^+)"^^ = 
^ (1 00) + 1 11)) are shared between Alice and Bob. Sup- 
pose Alice uses an n-fold Clifford encoder U to encode 
a /c-qubit state |?/^) in n physical qubits (including the c 
halves of the ebits on Alice's side) and then sends it to 
Bob. This is called an [[n, /c, d] c]] EAQEC code for some 
minimum distance d. The encoded state is 
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where the superscript A or B indicates that the operator 
acts on the qubits of Alice or Bob, respectively. Let g^ j = 
UZjU^ and h^j = UXjU^ for j = 1, • • • , n. The encoded 
state is stabilized by g\ (g) for i = A: + 1, • • • , /c + c, 
/I'f (g) Xf for i = A: + 1, • • • , A: + c, and g'f (g) for i = 
/c + c+1, • • • , n. These g'^^s and h'j^s are called simplified 
stabilizer generators. (We will omit the superscripts A, B 
as there is no ambiguity.) Note that g^i and h\ anti- 
commute with each other and they commute with all 
other generators g'j and h'j for j ^ i. Thus g'i and h'i 
are symplectic partners for i = 1, • • • , /c + c. An EAQEC 
code is defined by the simplified stabilizer group of the 
encoded state: 

— k+l-> ' ' ' •> 9' k+c-> k+1-, ' ' ' h' kJ^c-, g' kJ^c+l-> ' ' ' •> 9' n)- 

The symplectic subgroup of S' is Ss = (^'/c+i, " ' 1 g' k+c^ 



h'k+i^ ••• , h'k+c)-, and the isotropic subgroup of S' is 
Si = (^/c+c+i,, ••• , g'n)' = Ss X Si, so that 5' is 
generated by the generators of Ss and Sj. The minimum 
distance of the EAQEC code is the minimum weight of 
any element in JV{S^)\Si [13 . The decoding process of 
an EAQEC code is similar to that of a standard stabilizer 
code. 



III. DETERMINING THE SYNDROME 
REPRESENTATIVES FOR EAQEC CODES 

For a channel with low error rate, we would like to 
choose the set of syndrome representatives T to contain 
lower- weight error operators, since these errors are more 
likely to occur. A simple algorithm to define a set T is 
as follows: 

(a) Let T = and = 

(b) Find an error operator E e Qn\{T[JT') with the 
lowest weight, and compute its syndrome s. 

(c) If there is no Eg G T, set Eg ^ E d^ndT ^ T[j{Es}. 

(d) If there is an Es G T, set T ^ V U {E}. 

(e) If \T\ < 2^-^, go to (b). Else, output T. 

This algorithm finds a set T of minimum weight, that is, 
it minimizes the quantity 



In the usual paradigm of EAQEC codes, it is assumed 
that Bob's qubits suffer no error. However, this assump- 
tion might not be true in practice. Suppose that Alice 
uses a noisy channel JVa to communicate with Bob and 
Bob's ebits suffer from a storage error channel JVb- As- 
sume both A/a and A/^ are depolarizing channels. Let Pa 
and pi) be the depolarizing rate of A/a and A/b, respec- 
tively. Define 
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for = 0, • • • , n, and 
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(2) 



(3) 



for 



0, • • • , c. An error operator Ea^Eb of JVa'^J^b 



occurs with probability Qwa^wb^ where Wa = wt{EA) and 
Wb = wt{EB). 

To correct some errors on Bob's qubits, we have to de- 
sign a quantum code such that these errors are either syn- 
drome representatives, or degenerate errors of other cor- 
rectable errors. It is more complicated to determine the 
syndrome representatives in the case of EAQEC codes. 
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We will show that this criterion leads to high channel 



FIG. 1. (Color online) Contour plot of ^2^0— ^0^1- The shaded 
region is where ^2^0 — Qori > 0. 



since the error probabilities are different on Alice's and 
Bob's qubits. For an EAQEC code to correct some errors 
on Bob's qubits, we must sacrifice some of its ability to 
correct channel errors. For example, consider Bowen's 
[[3, 1,3; 2]] code with the following stabilizer generators: 



X z z 

Z Z X 

Z Y Y 

Y Y Z 



X I 

I X 

Z I 

I z. 



(4) 



Alice's qubits and Bob's qubits are on the left and the 
right of the vertical line, respectively. The error operators 
X4 and Y1Y2 have the same error syndrome. X4 is an 
error operator on Bob's side with probability 



qori 1 



1 



If Bob's qubits are error free, the code can correct the 
weight-2 error Y1Y2 on Alice's side, which occurs with 
probability 



fidelity in Section [Vl| We propose two coding schemes in 
the next two sections. 



IV. EAQEC CODES THAT ARE EQUIVALENT 
TO STANDARD STABILIZER CODES 

The stabilizer generators of the five qubit code [12 are: 

X Z Z X I 

I X Z Z X 

X I X z z 

Z X I X z 

After some row operations, these are equivalent to the 
stabilizer generators of Bowen's [[3, 1, 3; 2]]^^ EAQEC 
code in Q. Thus Bowen's quantum code is equivalent to 
the [[5,1,3]] standard stabilizer code, and has the same 
error-correcting ability on Bob's halves of the ebits as on 
Alice's qubits. We use the subscript AB to indicate an 
EAQEC code with this property. If we move the vertical 
line and put it between the fourth and fifth qubits in Q , 
we obtain the stabilizer generators of the [[4, 1, 3; 1]]ab 
EAQEC code in [13]. These observations inspire us to 
find EAQEC codes that are equivalent to standard sta- 
bilizer codes and the answer is straightforward as follows. 
Given a set of stabilizer generators {^1,^2, ' ^9n-k} of 
an [[n, /c,(i]] stabilizer code, after Gaussian elimination 
(and reordering qubits as necessary), they can be writ- 



g'c ^ ^c, 



/ for some c < [^^^^J such 



ten as g'^ 

K^Xc, ^/c+i^^, • • • , yn-k-c^-' ^^^^^^ - - L 2 

that the simplified generators g[ and /i- anti-commute 
with each other, and they commute with other simpli- 
fied generators ^ /ij for z = 1, • • • , c and j, / 7^ z; while 
g'c+i-, • • • , g'n-k-c commute with all the simplified genera- 
tors. Consequently, {g'l^ ■ ■ , g'^_j^_^^ h[^ - • ■ ^h'^} defines 
an [[n — c, /c, c]] EAQEC code. We summarize the above 
result in the following theorem and provide a lower bound 
on c. 

Theorem 2. Suppose H = [Hx\Hz] is the check matrix 
of an [[n,/c,(i]] standard stabilizer code. After Gaussian 
elimination, H can be written in the following standard 
form: 



q2ro = 1 



;Pa 



1 

-^Pa 



;Pb 



We can instead choose to correct X4 if qori > ^2^0- This 
is a tradeoff between correcting channel errors or storage 
errors. We plot q2ro — qori as a function of Pa and pb in 
Fig. [1] It can be seen that Y1Y2 is a more likely error 
than X4 when is small or Pa is large. 

To sum up, the set of syndrome representatives T 
should be chosen to minimize 

EeT 
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B ^sxs 
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F 
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where < s < n — k. Then there exists an [[n — 
c^k^d; c]]ab EAQEC code for some < c < 5 such that 
the code can correct J -qubit errors on either Alice's 
or Bob's qubits. In the case of nondegenerate quantum 
codes, s is bounded by 

k, 



d-l<s<l- 
Proof. An [[n — c, /c, c]]ab EAQEC code is defined by 
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C 


B Isxs 



the check matrix H with the last c columns of both Hx 
and Hz being removed. 

If d = 1, the result is trivial. Assume d > 2. We per- 
form Gaussian elimination on the check matrix H to pair 
the rows in the form ([5|. When the process of Gaussian 
elimination stops, either all the rows of H are paired or 
some rows are left unpaired. In the first case, 



H 



We have 2s = n — k > 2{d—l) by the quantum singleton 
bound, and thus s > d—1. In the second case, we can find 
a column in Hx (or Hz) that has a '1' at the (2s + l)-th 
position and O's elsewhere, and the corresponding column 
in Hz (or Hx) is of the form 

[uvw 0---0]^, 

where u and v are two binary s-tuples and i(; is or 1. 
Then H is in the following form: 



H 



1 



A' 



1 

CO 

" " " 1 


E' . 




B' 



1 



1 



w 


F' . 




Note that a subset of the last (s + 1) columns of Hx and 
the last (5 + 1) columns of Hz are linearly dependent. 
Thus we can find an element of weight at most (5 + 1) 
corresponding to these columns that is in the normalizer 
group of the stabilizer group. In the case of nondegen- 
erate quantum codes, the minimum distance is the min- 
imum weight of an element in the normalizer group of 
the stabilizer group. Since the minimum distance of the 
quantum code is (i, this implies that 



> d-l. 
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By this theorem, we can "move" some ancilla qubits 
to Bob's side for any nondegenerate stabilizer codes and 
obtain an EAQEC code. The case c = d — 1 of the 
above theorem is also observed by Wilde and Hsieh [15] 
from the viewpoint of purification and tracing qubits [28] . 
Note that Shor's [[9, 1, 3]] degenerate code [l] and Bacon's 
[[n^, l,n]] degenerate codes ^ also satisfy this theorem 
(by pairing up with Sf). It is conjectured that all 
degenerate codes satisfy the lower bound in Theorem [2] 
We have checked that all the optimal quantum codes in 



Grassl's table [30^ satisfy the lower bound in this theo- 
rem. 

EAQEC codes that are equivalent to standard stabi- 
lizer codes can correct errors on the qubits of both Alice 
and Bob. The decoder of the corresponding standard 
stabilizer code can be adopted to decode these EAQEC 
codes. These codes may perform better in practice than 
their corresponding standard stabilizer codes, for there 
are fewer physical qubits transmitted through the noisy 
channel, and the storage error rate is generally lower than 
the noisy channel error rate. 

In the case of CSS codes, the standard form of a parity 
check matrix is 



H 



^ ^sxs 







B L 





{n—k — s)x{n—k — s) 



Consequently, we have the following theorem. 
Theorem 3. An [[n, CSS code, obtained from an 



[n, /c', d] classical dual-containing code with k = 2k' 
gives [[n — c, /c, c\]ab EAQEC codes for < c < 
n-k'. 



n. 



This shows that any CSS codes can be transformed into 
EAQEC codes that correct errors on both Alice's and 
Bob's qubits. The decoding method in this scheme is 
exactly the same as that of standard CSS codes. Since 
many stabilizer codes are based on the CSS construction, 
we can take advantage of these codes while also having 
the power of entanglement. 

Example 1. From the [[7,1,3]] Steane code, we ob- 
tain a [[4, 1, 3; 3]]ab, a [[5, 1, 3; 2]]ab, and a [[6, 1, 3; 1]]ab 
EAQEC code. The [[4, 1, 3; 3]]ab and the [[6, 1, 3; 1]]ab 
EAQEC codes were also found by Wilde and Fattal [19| 
and Shaw et al [20 , respectively: 



XXIX 
XXXI 
X I X X 
Z Z I z 
Z Z Z I 
Z I z z 



XII 
I X I 
I I X 
Z I I 
I Z I 
I I z 



The quantum singleton bound says that for an [[n, /c, d]] 
quantum code, 

n-k>2{d-l). 

If the parameters n, /c, d of a standard stabilizer code 
achieve the quantum singleton bound, or n — k = 2((i— 1), 
then EAQEC codes equivalent to these standard stabi- 
lizer codes will achieve the singleton bound for EAQEC 
codes: 

n + c- /c < 2((i - 1), 

and we have the following theorem. 

Theorem 4. A standard stabilizer code that achieves 
the quantum singleton bound gives an [[n — c, /c, c\]ab 
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EAQEC code for some c that achieves the singleton 
bound for EAQEC codes. 

Example 2. The [[3, 1, 3; 2]]ab and [[4, 1, 3; 
EAQEC codes derived from the five qubit code achieve 
the singleton bound for EAQEC codes. 

Example 3. Using MAGMA [31^ to find optimal stan- 
dard stabilizer codes, we obtain several optimal EAQEC 
codes that achieve the linear programming bounds in [17] 
by Theorem [2] 

[[15, 10, 4; 5]]ab, [[14, 11, 3; 3]]ab, [[13, 9, 4; 4]]^^, 
[[13, 10, 3; 3]]ab, [[12, 9, 3; 3]]^^, [[11, 8, 3; 3]]^^, 
[[10, 6, 4; A\]ab. [[10, 7, 3; 3]]ab, [[9, 6, 3; 3]]ab, 
[[7, 4, 3; Z]]ab. [[8, 4, 4; 4]]^^, [[6, 2, 4; A]]ab. 
[[7,3,3;1]]ab,[[6,3,3;2]]^b,[[6,1,5;5]]ab, 
[[4,1,3;1]]ab,[[4,1,3;3]]ab,[[3,1,3;2]]^b. 

We also found the following EAQEC codes that have the 
highest minimum distance for fixed n and k to the best 
of our knowledge: 

[[15, 4, 8; 11]]ab, [[13, 5, 6; 8]]ab, [[12, 6, 5; 6]]ab. 



QUANTUM CODES WITH TWO 
ENCODERS 



In the previous section, we discussed EAQEC codes 
that are equivalent to standard stabilizer codes. Most 
optimal EAQEC codes are not equivalent to any stan- 
dard stabilizer codes, such as the entanglement-assisted 
repetition codes [HI [17]. We would like to exploit the 
high error-correcting ability of these quantum codes even 
in the presence of storage errors on Bob's side. This can 
be achieved by using another quantum code to protect 
Bob's qubits. 

Assume Alice uses the encoding operator Ua of an 
[[n, /c, (i^; c]] EAQEC code to protect her information 
qubits. Suppose also that there are m — c > ancilla 
qubits on Bob's side, and that Bob applies the encoding 
operator Ub of an [[m,c, d^]] standard stabilizer code to 
protect his c qubits. The encoding operator on the whole 
is Ua^Ub' We use the notation [[n, /c, d^; c]] + [[m, c, cIb]] 
to represent such a composite quantum code. If there are 
no ancillas on Bob's side, the set of stabilizer generators 
is equivalent to that of an EAQEC code for some encod- 
ing operator [/^ (g) after Gaussian elimination. 

The initial state is (g) |0)^^-^-^ (g) |0)^^-^ 
The encoded state has the following stabilizer generators: 
gi (g) ui, hi g) vi, g2 g) 1^2, h2 V2, - - - , 9c^ ^c, ^ ^c, 
Qc+i g) /, • • • , gn-k-c <S) I, I ^ • • • , I ^Um, whcrc 



UBZf'Uj, = Ui and UbX^U], 

A straightforward decoding process of the 
[[n, /c, dA] c]] + [[m, c, (Ib]] quantum code is that Bob first 
decodes his c ebits in the [[m,c, d^]] quantum code, 
and then he decodes the k information qubits hiding 



in the [[n, /c, cZ^; c]] EAQEC code. Or we can treat the 
combination code as an [[n + m^k^dc]] code, which has 
a more complicated decoding circuit but a potentially 
higher error-correcting ability. 

Under what condition is an [[n, /c, c]] EAQEC code 
not equivalent to an [[n + c^k^d]] standard stabilizer 
code? An [[n, /c, c]] EAQEC code satisfies the Ham- 
ming bound for EAQEC codes [21]: 



< 2 



n— fc+c 



for t = [^^J. If the parameters do not satisfy the 
quantum Hamming bound for (nondegenerate) quantum 
codes: 



t 



3^ 



J 



-k+c 



an [[n + c, /c, d]] code does not exist. 

Example 4. Consider the [[7, 1,5; 2]] code in [M]. We 
have E5=o3^(i) = < 2^-^+^ = 256. However, 

E5=o3^ rr) " ^ 2^+^-^ = 256. Hence there is 
no [[9, 1, 5]] code. 

However, the Hamming bound is not tight. For ex- 
ample, the [[7, 1, 5; 3]] code is not equivalent to any stan- 
dard code, but the parameters [[10, 1,5]] satisfy the quan- 
tum Hamming bound. A better bound, such as the lin- 
ear programming bound, can be applied here. If the 
parameters n + c, /c, d violate any upper bound on the 
minimum distance of the quantum code, then such an 
[[n + c, /c, d]] standard stabilizer code does not exist, and 
the [[n, /c, (i; c]] EAQEC code is not equivalent to any 
standard stabilizer code. We can check the result for 
n + c < 30 from the tables of stabilizer codes in [TT] and 

Several EAQEC codes that not equivalent to standard 
codes were found in [14]: [[n, l,n;n — 1]] for n odd, 
[[n,l,n-l;n-l]] forn even, [[7, 1,5; 2]], [[7, 1,5; 3]], 
[[7, 2, 5; 5]], [[9, 1, 7; 4]], [[9, 1, 7; 5]], [[9, 1, 7; 6]], [[9, 1, 7; 7]], 
[[8, 2, 5; 4]], [[8, 3, 5; 5]], [[13, 3, 9; 10]], [[13, 1, 11; 11]], 
[[13, 1,11; 10]], [[13, 1,9; 9]], [[13, 1, 9; 8]]. Also, the 
[[15, 7, 6; 8]], [[15, 8, 6; 7]], [[15, 9, 5; 6]] EAQEC codes are 
not equivalent to any known standard quantum stabilizer 
code. 

Example 5. The [[n, l,n;n - 1]] EAQEC codes for n 
odd saturate the quantum Hamming bound with equality 
as follows. The number of correctable X 01 Z errors is 
l^^^J number of correctable error syndromes is 



22(n-l) 



Suppose Alice uses the [[5, 1,5; 4]] entanglement-assisted 
repetition code and Bob applies the optimal [[10,4,3]] 
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X 



FIG. 2. Illustration of the weight distributions of the syn- 
drome representatives of the two decoding method. 



quantum code to protect his 4 ebits. Then the whole 
quantum code [[5, 1,5; 4]] + [[10,4,3]] can protect two 
channel noise errors on the 5 qubits that Alice sends 
through the channel, and one storage error on Bob's 10. 
On the other hand, the optimal quantum code using 15 
qubits to encode one information qubit is the [[15,1,5]] 
quantum code ^], and it can correct an arbitrary two 
qubit error. 

Compared to an [[n + m, /c, d]] standard stabilizer code, 
the number of qubits going through the noisy channel is 
much less using the [[n, /c, d^; c]] + [[m, c, (i^]] quantum 
code. Consequently, as long as rate the storage error 
is reasonably small compared to the channel error rate, 
the [[n, /c, (i^; c]] + [[m,c, d^]] quantum code has better 
efficiency, while keeping the same error-correcting ability 
against the channel noise. The performances of different 
coding schemes will be compared in Section |VII| 

For combination EAQEC codes, the syndrome repre- 
sentatives can be chosen as in the case of standard sta- 
bilizer codes if we treat the code as a single stabilizer 
code. If we use two sequential decoders, we choose two 
sets of syndrome representatives Ta, Tb as in the case 
of standard stabilizer codes. Observe that T = Ta x Tb 
is the set of syndrome representatives of the combination 
code. The weight distributions of the syndrome represen- 
tatives of the two decoding method are illustrated in Fig. 
[2] The X-axis and y-axis represent the weights on Bob's 
and Alice's qubits, respectively. The weight distribution 
of the syndrome representatives using two sequential de- 
coders is always a rectangle. If the error probabilities 
are the same on Alice's and Bob's qubits, the weight dis- 
tribution of the syndrome representatives using a single 
decoder looks like a triangle. If the error probabilities are 
different on Alice's and Bob's qubits, the shape varies ac- 
cording to the error probabilities. Note that the area of 
the triangle is equal to the area of the rectangle. 

VI. CHANNEL FIDELITY 

A. Formula for the Channel Fidelity over the 
Depolarizing Channel 

Let f be a quantum channel operating on the input 
state p, which lies in a state space H' of dimension m. 



and the output state S{p) also lies in . Suppose the 
quantum channel £ has the operator-sum representation 
£{p) = ^^EipE\^ where the operation elements {Ei} 
satisfy T^iElEi = I- Let |V) = where 
{\j)} is a basis of the input space H' . The channel fidelity 
[241 of f is defined as 

Fc{£) ={^\ {1 £\^)m m = ^T. 

i 

The channel fidelity can be used as a measure of the 
performance of a quantum error-correcting code over a 
noisy channel. Suppose the /c-qubit information state 
I?/;), after being encoded by the encoding channel U : 
y^^k _^ ^ is transmitted through a noisy channel 
M : n®"" H^^, and then decoded by the decoding 
channel V : ^ U®^ . Then FciVMU) serves as 

a merit function of the quantum code with encoding- 
decoding pair {U^ V} over the noisy channel JV. 

The encoding channel U can be written as 

U{\^){^\) = /7£;(/^^0|O)^"-^)|V^)(V^|(/^^0(O|^"-^)/7];, 

for some unitary Clifford encoder Ue- We assume there 
are (n — k) ancillas and we implicitly use the equation 

for any /c-qubit state |?/^). The decoding channel V con- 
sists of the following steps: syndrome measurement {Pg}, 
correction {Cg}, decoding U^, and partial trace of the 
ancillas, which have been explained in Section |n[ The 
overall process of the decoding channel is 

I s 

where {|/)} is the standard basis for the state space of 
the (n — k) ancilla qubits, and |/) has a binary repre- 
sentation 1/1/2 •• • In-k)' Thus V has operation elements 
{(/^^ (g) {l\)UljCsPs}i,s' Suppose the noisy channel J\f 
is the independent n-fold depolarizing channel with 
operation elements {s/piEi}^ where Pi is defined in ([T]). 

The composite channel VMU has operation elements 

J l,i 

and its channel fidelity is 

l,s,i 

There are a total of 2^-^ • 2^-^ • 4^ = 4^^-^ terms in the 
sum, and each term is a product of one (2^ X 2^) matrix. 
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five (2^ X 2^) matrices, and one (2^ x 2^) matrix. Thus 
the complexity of the calculation of the channel fidelity is 
Q{n ' 43^-2^^ complexity of multiplication of (2"^ x 
2^) matrices is r^(n2^)), which is almost impossible to 
calculate for n > 10. We will show how to reduce the 
complexity to something more manageable. 

Lemma 5. If the error syndrome of Ei is 5^, then 

|trW,i,.)|' = 

for s Si. In addition, the measurement result is Si with 
certainty. 



Type of codes 


Channel fidelity 


bit-fiip code 


2^ 


[[4,2,2]] code 


1- ip+ iP" 


[[5,1,3]] code 

LL ' ' JJ 




[[7,1,3]] code 


1 147 2 ,189 3 1785 4, 1155 5 


[[8,3,3]] code 


1 245 2 1 1449 3 8U29 4 , 12743 5 

2961 J , 763 21 8 
128 ^ ^ 128^ 32^ 


[[9,1,3]] code 


1 - V + - + 

567^6 1 225^7 27^8 , 154S 9 
32 ^ 32 ^ 16^ ^ 8192^ 



TABLE I. channel fidelity of different quantum codes in the 
depolarizing channel 



It is straightforward to check that PgEiUEi 
|Q^(8)n-/cj ^ Q ^Yie facts that Pg's are orthogonal 

to each other and Pq is the projector on the code space, 
and the above lemma follows naturally. Next we show 
that only the error operators in T x 5 have nonzero con- 
tribution to the channel fidelity. 



Lemma 6. 



1 

2^ 



if G T X 5; 
otherwise. 



Proof. It is straightforward to verify the lemma for the 
case that Ei G T x S. Now assume Ei ^ T x S and we 
have Cs^Ei G JV{S)\S. Since t/^; is a Clifford unitary op- 
erator. Ml 0M2 — UljCsiEiUE G 0n is a Pauli operator, 
where Mi G Qk and M2 G Qn-k- Since Ei ^ T x we 
have Ml not equal to the identity and tr(Mi) = 0. Let 
the matrix representations of Mi and M2 be [aij] and 
[bij]^ respectively. Then 



^ |tr ((/^^ {1\)uICs,Ps,E,Ue{I^^ ^ |0)®"-^)) 
I 

= ^ |tr ((/^^ {l\)Mi M2(/^^ |0)^^-^)) 1' 
I 

= ^|5M-tr(Mi)|'=0, 



where the second equality follows by explicitly writ- 
ing down the matrix multiplications of {I^^ (/|)Mi 

□ 



Proof. 



Fc(I>A/-Zi)=^^|trW,,,..)r 

14 

= ^ 

EieTxS 
n 

^ ^ ^wQw 1 
w=0 



(6) 



where is the number of G T x 5 of weight k;, and 
is defined in (!]). The first equality follows by Lemma 



|5]and the second equality follows by Lemma [6] 



□ 



There are 4^~^ terms in ^ and each term is gener- 
ated by vector addition. The complexity is now reduced 
to r^(n-4^~^). From the above theorem, we find that the 
channel fidelity for the depolarizing channel is the proba- 
bility of correctable errors. For a single information-qubit 
code, the channel fidelity is the probability that the in- 
formation qubit can be correctly recovered. The formula 
for the channel fidelity of several quantum codes is shown 
in Table. [D 

The channel fidelity of the quantum code defined by a 
stabilizer group S over a depolarizing channel depends on 
the syndrome representatives T and the stabilizer group 
5, and it can be optimized over the choices of T. (This 
is to optimize over different decoding schemes.) 



B. Channel Fidelity for EAQEC Codes 



We derive the following theorem with the help of 
Lemma [5] and Lemma [6l 

Theorem 7. The channel fidelity of a quantum code, 
with a stabilizer group S and a set of syndrome repre- 
sentatives T, over the depolarizing channel, is the weight 
enumerator of the probability distribution {qw} of the 
elements in T x 5. 



Suppose that Alice uses a noisy channel Ma to com- 
municate with Bob and Bob's ebits suffer from a storage 
error channel Mb- 

Definition 8. The channel fidelity [231 [24] of an EAQEC 
code with encoding and decoding processes {U^ V} over 
the noisy channel Ma ^ Mb is 

Fc(V{Ma®Mb)U). 
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. 10-= 

FIG. 3. Comparison of two [[3, 1,3; 2]] EAQEC codes in 
terms of channel fidelity. The [[3, 1, 3; 2]]ab code performs 
better in Region B, while the repetition code performs bet- 
ter in Region A. The region for < 0.015 and pb < 0.2 is 
enlarged on the right. 



Suppose an [[n, /c, c]] EAQEC code is defined 
by a simplified stabilizer group = {g[r ' ' ^9n-k^ 
^ir*' 7^c) ^ Sn- We extend the simplified stabi- 
lizer group to a stabilizer group S = (^i,--- ^gn-k^ 
hi," - , he) ^ Qn+c, where Qi = g'^ ® Zi, hi = h[ X^, 
for i = 1, • • • , c, and = (g) / for j = c + 1, • • • ^n — k. 
We choose a set T of syndrome representatives for the 
noisy channel Ma ^Mb- Let Pa and p5 be the depolariz- 
ing rate of Ma and Mb , respectively. Now we can apply 
Theorem [7| to find a formula for the channel fidelity of 
EAQEC codes. 

Theorem 9. The channel fidelity of an EAQEC code 
with a stabilizer group S and the set of syndrome rep- 
resentatives T over the depolarizing channel Ma ^ Mb 
is the weight enumerator of the probability distribution 
{Qw'^w'} of the elements in T x 5, where p^ and r^f are 
defined in ([2| and (|3|. 

Note that the channel fidelity can be optimized over 
the choice of T. Since the two noisy channels Ma and Mb 
have different error rates, error operators on these two 
channels should be differently weighted. In the extreme 
case of p5 = and ro = 1, the error-correction condition 
for EAQEC codes says that {Ei} is a set of correctable 
error operators if E^Ej ^ M{S^)\Si, where Si is the 
isotropic subgroup generated by {g^-^-i^ ' • • ^g^-k}- Note 
that ro = 1, and = for z 7^ 0. We can similarly 
formulate the channel fidelity. 

Corollary 10. The channel fidelity of an EAQEC code 
with a stabilizer group S and a set of syndrome represen- 
tatives T over the depolarizing channel Ma ^ Ib is the 
weight enumerator of T x 5/, which is a polynomial in 

{Qw}- 

Therefore, we would like to choose T to consist of likely 
error operators from the noisy channel Ma- 

Example 6. We compare the channel fidelities of the 
[[3, 1,3; 2]] entanglement-assisted (EA) repetition code 



0.9 - 




FIG. 4. (Color online) Contour plot of the difference between 
the two fidelities: Fc{VW {T®J' ®T®'')Ur) - Fc{Vi'' {TZ"" ® 
^®2-^^ABy ^jjg shaded region is where Fc{W{Tp^J ® 



Ha, Fc{V-s%%f®Tpf)W), and Bowen's [[3, 1, 3; 2]]ab 
quantum code [21_, FciV^^ {T®^ ®T®'^)Ui^). Bowen's 
[[3, 1, 3; 2\\ab quantum code is equivalent to the five 
qubit code and has T = {I, Xi,--- jX^, Zi,--- ,^5, 
Yi,''' ,^5}, while the [[3, 1,3; 2]] EA repetition code is 
designed under the assumption that Bob's qubits are per- 
fect and hasT = {/,Xi,-- - ,X3,Zi,--- ,Z3,Yi,--- ,¥3, 
XiZ2,XiZ3, ZiX2,X2Zs, ^1X3,^2X3}. We have 

Fc{Vi^{Tpf ® rf,')U^'') = qoro + Qqiro + Gqsro + 6qon 
+ 36q2ri + BAq^ri + I8qir2 + Slq2r2 + 4593r2, 

and 

Fc{Vl%Tpf TpfMl = Qoro + 9q,ro + 6q2ro 
+ ISqiTi + 38(72^1 + 40(73ri + 18qir2 + 55(72^2 + 71q3r2. 

The channel fidelities of these two EAQEC codes are 
compared in Fig. |3] and Fig. [4j In Fig. [3j the curve 
of the boundary between the two regions passes the ori- 
gin. The region A in Fig. [3] corresponds to the shaded 
part in Fig. [4] In region B, Bowen's code is better than 
the EA repetition code. 

In the extreme case of = 0, we have 

Fc(pr(7;f ® I'^^M'') = 90 + 9</i + 693, 

and 

FciWiTZ'' ® I^^Ml =qo + 9qi + 692- 

The EA repetition code corrects more lower-weight er- 
rors, and hence it has higher channel fidelity. 

The channel fidelity of the [[n, /c, d^; c]] + [[m, c, d^]] 
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quantum code depends on the decoding process. If we 
treat the combination code as an [[n + m^k^dc]] code, 
its channel fidelity can be computed as Theorem |9] If 
we use two sequential decoders, it's different. First, the 
4^ possible error operators on Bob's side collapse to 4^ 
logical errors. Each can be obtained from error 
operators, after the decoding process of the [[m,c, (i^]] 
quantum code. If the 4^ errors occur uniformly, the 
decoded errors also occur uniformly. However, this is not 
the case for the depolarizing channel Afs- the 4^ errors 
on Bob's ebits occur according to a distribution {r[^,} 
that depends on the decoding process. We can find the 
channel fidelity as in Theorem [9j except that the errors 
on Bob's ebits follow a new distribution {Qw^w'}' 
will compare the channel fidelities of these two decoding 
methods in Section IVlII 

C. Approximation of Channel Fidelity 

The number of terms in the formula for channel fidelity 
over the depolarizing channel is 4"^+^-^^ which grows ex- 
ponentially in n + c — /c. It is difficult even to build a 
look-up table for decoding when n is large. However, it 
is possible to approximate the channel fidelity efficiently. 

The channel fidelity of a quantum code over a depolar- 
izing channel is the probability that the decoder output 
is correct, which can be lower bounded by 

Pr( {correctable errors}) 

> Pr({syndrome representatives}), 

or 

Pr({correctable errors}) 

> Pr( {errors of weight less than or equal to L^T^J})^ 

where d is the minimum distance of the quantum code. 
When the depolarizing rate is low (< 0.2), these bounds 
are fairly tight and are good approximations of the true 
channel fidelity. Dong et al. defined an "infidelity" func- 
tion to characterize the performance of quantum codes 
[3? , which is an approximation of the channel fidelity in 
the case of k = 1. 

We can also apply Monte Carlo methods to approxi- 
mate the channel fidelity [25 , especially when the num- 
ber of physical qubits involved is large. The steps of the 
simulations are as follows: 

1. Fix a depolarizing rate p. 

2. Randomly generate an error operator E accord- 
ing to the probability distribution of a depolarizing 
channel. 

3. Compute the error syndrome and apply the correc- 
tion operator and decoding operator. 

4. If there is no logical error after decoding, E is cor- 
rectable. 




0.2 - 



0.1 - - 

qI < < < < 

0.2 0.4 0.6 0.8 1 

Depolarizing rate p 

FIG. 5. (Color online) The approximations of the channel 
fidelity of the [[11,1, 5]] quantum code. 
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lower bound of [[24,1,8]] 
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FIG. 6. (Color online) The approximations of the channel 
fidelity of the [[24, 1, 8]] quantum code. 



5. Repeat steps 2 to 4 times. 

6. Output the channel fidelity as the number of cor- 
rectable errors in the experiment divided by N. 

Two applications of the Monte Carlo method to the 
[[11, 1, 5]] and [[24, 1,8]] codes obtained by MAGMA [31] 
are shown in Fig. |5]and Fig. [6j In Fig. [Sj the exact chan- 
nel fidelity of the [[11,1,5]] code is also plotted, and it 
can be observed that the simulations quickly converge to 
the exact channel fidelity. On the other hand, in Fig. [6j 
only a lower bound of the channel fidelity of the [[24, 1, 8]] 
code is given, since computation of the channel fidelity 
is difficult (4^^ terms in the formula). However, Monte 
Carlo simulations also converge quickly from A/" = 10^ to 
10^ points, which are much less than 4^^ ^ 7 x 10^^. 
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VII. PERFORMANCE ANALYSIS 

In this section, we compare the channel fidehty of the 
[[5, 1,5; 4]] + [[10,4,3]] quantum code with other quan- 
tum codes. A comparable standard stabilizer code with 
n = 15 and /c = 1 is a [[15,1,5]] stabilizer code, while 
the smallest quantum code with d = 5 is a [[11,1,5]] 
stabilizer code [30 . By Theorem [2j we can obtain a 
[[6,1,5;5]]ab EAQEC code from the [[11,1,5]]. The 
[[5, 1,5; 4]] EAQEC code is also shown as a reference. 
All the quantum codes above encode k = 1 information 
qubit and hence can be compared with no ambiguity. 
The channel fidelity of these quantum codes are plotted 
in Fig. [7| For simplicity, the [[5, 1, 5; 4]] + [[10, 4, 3]] quan- 
tum code using a single decoder is treated as a standard 
stabilizer code, which means that the weight distribution 
of the syndrome representatives is like a triangle as in 
Fig.|2j 

When pb = OMpa, the [[6, 1, 5; 5]]ab EAQEC code 
is the best code for Pa < 0.38, and [[5, 1,5; 4]] + 
[[10,4,3]] quantum code with sequential decoders and 
the [[5, 1, 5; 4]] EAQEC code have better performance for 
Pa > 0.38. The [[5, 1,5; 4]] + [[10,4,3]] quantum code 
with a single decoder has a comparable fidelity to the 
[[6,1,5;5]]ab EAQEC code. 

As the ratio of p^ to Pa decreases, the perfor- 
mance of the [[5, 1,5; 4]] + [[10,4,3]] quantum code and 
the [[6, 1,5;5]]ab, [[5, 1,5; 4]] EAQEC codes get better. 
When pk = 0.5pa, the [[5, 1,5; 4]] + [[10,4,3]] quantum 
code, using either decoding method, performs better than 
the [[6,1,5;5]]ab EAQEC code and the [[11,1,5]] stabi- 
lizer code. When the ratio of pb to Pa is below 0.35 in the 
last two cases, the [[5, 1, 5; 4]] + [[10, 4, 3]] quantum code is 
the best choice. The [[5, 1,5; 4]] EAQEC code beats the 
[[11,1,5]] stabilizer code when = O.Olpa as expected. 
From these plots, the quantum codes that have the best 
performance are the [[6, 1, 5; 5]]ab EAQEC code or the 
[[5, 1, 5; 4]] + [[10, 4, 3]] quantum code. 

Fig. [sjplots the difference {Fsingie — Fseq) between the 
fidelity of a single decoder {Fsingie) and the fidelity of two 
sequential decoders (Fgeq) of the [[5, 1,5; 4]] + [[10,4,3]] 
quantum code. For low noise rates sequential decoding 
is not as good; however, it is easier to be implemented. 
Note that, in the simulations, the syndrome representa- 
tives using a single decoder are chosen for the same error 
probabilities on Alice's and Bob's qubits. The channel 
fidelity of the combination using a single decoder will be 
better if we optimize it over the choices of the syndrome 
representatives according to the different error probabil- 
ities. 



VIII. EAQEC AND ENTANGLEMENT 
DISTILLATION 

Up to now, we have assumed that Alice and Bob have 
perfect ebits before the communication, and any noise 
results from storage errors. In a more general situation. 



the exchange of perfect maximally-entangled states might 
not be possible, and entanglement distillation is needed. 
Our next direction is to find good strategies against this 
problem. 

We first introduce the entanglement distillation proto- 
col. Suppose Alice has the ability to prepare n pairs of 
maximally-entangled states (Bell states) 

1*+)^'" = (y^(ioo) + 111)) j = TP E l^)^!^)""' 

where i is the binary representation of the numbers be- 
tween and 2"^ - 1. The state |<|)^)^^ has the following 
property. 

Lemma 11. For any operator M on an m-qubit maxi- 
mally entangled state, 

/^|<|>+)^^ = (M^)^|<l>+)^^, 

where is the transpose of M. 

Suppose Alice uses an [[m, c]] stabilizer code defined by 
a stabilizer group S with generators /i, • • • , /m-c, and T 
is a set of syndrome representatives corresponding to S. 
Let /7 be a Clifford encoder of the stabilizer code, and 
\j) L for j = 0, • • • , 2^ — 1 be the logical states. The 
encoded |<l>+)^^ is 

j=o 

We know that {i\LEs-,Es^\j)L = Ss-,,sji,j and {Es\i)L} 
is a set of orthonormal basis vectors of 7/®". In the case 
of UU^ = I, we have 

^ EseT \j=o I 

If UU^ 7^ /, Alice applies the operator UU^ on half of 
the ebits. From Wilde's encoding circuit algorithm [33], 
an encoding operator can be implemented by a series of 
Hadamard gates, CNOT gates, SWAP gates, and phase 
gates. If phase gates are not used in the circuit, the 
circuit will satisfy IJU^ = I. 

Alice sends half of the ebits to Bob through a noisy 
channel Mc with depolarizing rate equal to Pc- The cor- 
rupted state is 

i?f|$+)«™ = ^'^' EfEfEf\j)i\j)l 

After performing a syndrome measurement, Alice obtains 
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FIG. 7. (Color online) Channel fidelity of different quantum codes for pb = 0.99, 0.65, 0.5, 0.35, 0.1, 0.01 pa- 



a syndrome a, which is a binary (m — c)-tuple, and she 
sends a to Bob through a noiseless classical channel. Now 
the state is 

where Ea is the syndrome representative of a. Let s{i) 
be the error syndrome of Ei. Bob measures the stabilizer 
generators /i, • • • , fm-c and obtains the syndrome 

b = a -\- s{i). 



The error syndrome s{i) can be retrieved by s{i) = a-\-b. 
He applies the correction operator and obtains the 
state 

E^Ef^,^EfE^\^+)r. 

Finally, they restore the state to the standard encoded 
state by applying the operator Ea^Ea^ followed 

by the decoding circuit U"^ U"^ and obtain the state 

if Ei is correctable. This is called an [[m, c]] entanglement 
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-[[10,2,7;8]] EAQEC code 
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FIG. 8. (Color online) Contour plot of the difference (Fgingie — 
Fseq) between decoding methods for the [[5, 1, 5; 4]] + [[10, 4, 3]] 
quantum code. The shaded region is where Fsingie— Fseq > 0. 



FIG. 9. (Color online) Comparison of the [[8, 3, 3]] distillation 
protocol plus [[5, 2, 3; 3]] EAQEC code, and the [[10, 2, 7; 8]]ea 
EAQEC code without distillation. The performance of the 
[[10, 2, 7; 8]]ab EAQEC code is better for pa < 0.45. 



distillation protocol. Let E denote the combination of the 
above processes. The entanglement fidelity [23 of this 
protocol is 



F(|<l>+)^^f ) = (|$+)^^(<l> 



The entanglement fidelity of the [[m, c]] entanglement dis- 
tillation protocol when the channel Afc is the depolariz- 
ing channel is just the channel fidelity of the correspond- 
ing [[m,c]] stabilizer code over the depolarizing channel 

After performing an [[m,c]] entanglement distillation 
protocol that produces c ebits with some entanglement 
fidelity, Alice can use an [[n, /c, c]] EAQEC code to send 
k information qubits to Bob. Assume Pb = and Pc = Pa 
for simplicity. The channel fidelity from the combination 
of an entanglement distillation protocol and an EAQEC 
code is similar to that of a combination code using two 
sequential decoders. On the other hand, Alice could di- 
rectly use an [[n^, /c, d'; e]] EAQEC code with the imper- 
fect ebits to send k information qubits to Bob without 
entanglement distillation. The channel fidelity of this 
process is just the channel fidelity of the EAQEC code. 

Example 7. Alice and Bob use an [[8,3,3]] entangle- 
ment distillation protocol to produce 3 perfect ebits from 
8 noisy ebits. Then Alice can use an [[5, 2, 3; 3]] EAQEC 
code to send quantum information to Bob. Or she can 
instead use a [[10, 2, 7; 8]] EAQEC code with the 8 cor- 
rupted ebits. Comparison of the channel fidelity of these 
two schemes is shown in Fig. [9] 

Example 8. Alice and Bob use a [[5, 1, 3]] entanglement 
distillation protocol to produce 1 perfect ebit from 5 noisy 
ebits. Then Alice can use a [[4, 1, 3; 1]] EAQEC code to 
send quantum information to Bob. Or she can instead 
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FIG. 10. (Color online) Comparison of the [[5,1,3]] dis- 
tillation protocol plus [[4, 1,3;1]] EAQEC code, and the 
[[6, 1, 5; 5]]£;a EAQEC code without distillation. The perfor- 
mance of the [[6, 1, 5; 5]] AS EAQEC code is slightly better for 

Pa < 0.11. 



use a [[6, 1, 5; 5]]ab EAQEC code with the 5 noisy ebits. 
Comparison of the channel fidelity of these two schemes 
is shown in Fig. [lO| 



In Fig. 11 



we plot the channel fidelity of four distil- 
lation protocols with the same [[4, 1, 3; 1]] EAQEC code. 
The channel fidelities of the 4 schemes with distillation 
protocols are better than the one without distillation 
protocol for p < 0.1, but the difference is small. The 
[[4, 1,3;1]] EAQEC code without distillation dominates 
the performance for higher p. 

In Fig. [121 we plot the channel fidelity of several coding 
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FIG. 11. (Color online) Comparison of different distillation 
protocols with the same [[4, 1, 3; 1]] EAQEC code. 
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FIG. 12. (Color online) Comparison of different combinations 
of distillation protocols and EAQEC codes that encode one 
information qubit. 



IX. DISCUSSION 

In this paper, we proposed two coding schemes for 
EAQEC codes when the ebits of the receiver suffer errors. 
In the first case we assume the ebits suffer storage errors. 
EAQEC codes that are equivalent to standard stabilizer 
codes have better performance than their corresponding 
stabilizer codes. Several such EAQEC codes are found to 
achieve the linear programming bound, and hence are op- 
timal. However, as long as the storage error rate is small 
enough, a quantum code with two encoders performs well 
if we start with an EAQEC code that is not equivalent to 
any standard stabilizer codes. We may choose the best 
quantum code according to the noise channel rate and 
the storage error rate in real situations. 

Any (nondegenerate) standard stabilizer code can be 
transformed into an EAQEC code by Theorem |2j Fam- 
ilies of quantum codes, such as quantum Reed-Muller 
codes quantum BCH codes [35 , 36 , quantum cyclic 
codes [371 EH], can be transformed into families of 
EAQEC codes. It is possible to construct EAQEC codes 
with a large number of information qubits but a small 
number of ebits that outperform standard codes. 

We developed a formula for the channel fidelity over 
the depolarizing channel, and used it to evaluate the per- 
formance of a variety of quantum codes. For large codes, 
the channel fidelity cannot be calculated exactly, but can 
be lower-bounded or approximated by Monte Carlo sim- 
ulations. A similar formula for the channel fidelity can 
be developed for other channels with only Pauli errors. 

We also compared EAQEC codes combined with en- 
tanglement distillation protocols to EAQEC codes de- 
signed to tolerate noisy ebits. It seems that EAQEC 
codes that can correct errors on both the qubits of sender 
and receiver can have better performance than the codes 
combined with an entanglement distillation protocol, at 
least for modest noise rates. For particular combinations 
of error rates and applications, it should be possible to 
optimize the choice of code to maximize the fidelity. This 
optimization is the subject of ongoing research. 



schemes that encode k = 1 information qubit. It can 
be observed that the EAQEC codes without distillation 
protocols have better performance for p < 0.2. 
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